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Abstract 

We study the free-propagation features of an optical field endowed with 
a non-uniform polarization pattern with elliptical symmetry. The fields 
derived in this way are called Elliptically Symmetric Polarized Beams 
(ESPB for short). Some properties of those fields are analysed. Moreover, 
it is shown how it is possible to obtain such light beams by applying the 
result to Bessel-Gauss beams. 



1 Introduction 

Recently there was a growing interest, in the Optics community, towards those 
optical fields possessing unusual polarization features; this because some phe- 
nomena have shown to be strongly dependent by vectorial properties of the field. 
For instance, it was experimentally demonstrated that a radially polarized field 
can be focused to spot size of extension much smaller than that obtainable by 
using a linear polarized one P^. Also, azimuthally and radially polarized fields 
[2] have found use in many sector of Optics, like particle beam trapping, optical 
microscopy, superrisolution pj], specially when those polarization features are 
associated to a particular intensity profile (the so-called donut beams are an 
example of that) which implies the presence of a strong longitudinal component 
of the optical electric field; moreover it was shown that beams possessing field 
lines in the shape of a spiral, namely spirally polarized beams, are suitable to 
describe fields with vortices in their transverse phase profile @]. In the present 
paper we wish to widen this scenario, by introducing coherent beams with a 
non-uniform polarization state with an elliptical symmetry. By doing so, we 
wish to show how some of aforesaid properties are easily extendable to other 
polarization's geometry and to emphasize the differences among all those polar- 
izations. The paper is organized as follows: in section 2, after having recalled 
briefly some definitions about cylindrical coordinate system, we specialize the 
subject to elliptical symmetry; in section 3 we derive some results in paraxial 
optics regime by applying the theory to Bessel-Gauss beams. 



2 Field and polarization symmetries 

Let us suppose we have a free-space beam, which propagates along a mean 
direction, say z-axis, on starting from a source-plane on which we know the field 
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distribution of amplitude and phase. In this context the term non-uniform, when 
addressed to the field polarization features, implicity refers to the transverse 
component of the field analysed through a rectangular (Cartesian) coordinate 
system. So we say that a polarization is non-uniform if, on moving along a 
fixed direction in the (x, y) plane, we see the vector field changing in direction 
when compared with the unit vector basis x and y. In reality, it should be still 
possible to find a new coordinate system respect to which the field mantains 
its polarization state. In other words, the field should be uniformly polarized 
respect to this new coordinate system. In a generical cylindrical orthogonal 
reference frame we can write the field, across the plane z = 0, as 

E = E T + E z z (1) 

where 

Er(gi,g 2 ,0) = E 1 (q 1 ,q 2 ,0)q 1 + E 2 ( qil q 2 ,0)q 2 (2) 

indicates the trasverse field and qi , q 2 and z indicate the vector basis of a gener- 
ical cylindrical orthogonal coordinate system, which is related to a Cartesian 
orthogonal system by the following relationships 

X = 9\{q\,q2) 

y = 52(91,92) (3) 

z = z 

where g\ and g 2 are smooth functions of q\, q 2 variables. Rectangular and 
circular ones are special cases of equations Q. In particular, speaking about 
circular symmetry, recently several authors have focused their attention on fields 
polarized in a radially or azimuthally way. In almost all these works also the field 
(i.e the amplitude and phase of each scalar component of the optical field, and 
not only the polarization) was choosen to be circularly simmetric. While this 
choice seems to be the most obvious and natural, anyhow it can be interesting 
to study the case in which this coupling in symmetry is removed, a case that 
does not have received attention at all. This fact is a little odd, specially if one 
thinks to the way in which a radially or azimuthally circular polarized field can 
be obtained in practice For instance, one of the methods to obtain a radially 
polarized beam is to combine coherently an Hermite-Gauss beam of order (0, 1) 
and an analogous beam, but of order (1,0), both arising from a laser beam, on 
supposing that amplitude and phase of the two beams are the same. There exist 
other metods, but all of them also suppose to have a very good control of the 
amplitude of the two component that one adds togheter. Actually this is true 
only approximately and in practice one can consider the effect of some difference 
in amplitude between the two emerging beams, a difference which gives rise just 
to an elliptical symmetry. In the next section we start by analysing the form 
of a generical field endowed with features that break both the circular than 
rectangular symmetry and we will discuss the relative properties. 

2.1 Vectorial beam with elliptical symmetry 

Let us suppose we have a field like that in equation and Ex be 

E T = Af(r, z) cos(9 + a)i + Bf(r, z) sin(0 + a)y (4) 
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where A and B indicate two arbitrary real constants, r and 9 are, respectively, 
the radial and angular coordinates of a circular-symmetric coordinate system, 
f(r, z) is a radial-circular function across a z plane and a is a constant angle. 
We suppose here that 4 / B. It is straightforward to show that the field in 
equation Q can be rewritten as follows 

Et = C/(r, z) [sinh £ cos(0 + a)x + cosh £ sin(0 + a)y] (5) 

by choosing the terms C and £ as 

C 2 = B 2 - A 2 (6) 

cosh£ = | (7) 

The first observation we do is trivial. If we fix the radial coordinate and we 
let to vary only 9- variable, we will see the arrow describing the electric field in 
equation (@J to describe an ellipse with semiaxis Af(r,z) and Bf(r,z). More 
interesting is to find the shape of the field lines for the aforesaid field. To achieve 
this goal we prefer to refer to a circular coordinate system. To obtain the field 
lines one needs to solve the equation 



dy E y Bsm{9 + a) 
dx E x Acos(9 + a) 

which becomes, in polar coodinates, 

sin Odr + r cos Odd sin 9 cos a + cos 9 sin a 

: = a : ; 

cos Odr — r sin 9d9 cos 9 cos a — sin 9 sin a 



(8) 



(9) 



where we have defined a new constant a = B/A. Equation can be trans- 
formed into the following relation 

dr [— (cos 2 9 + asm 2 9) cosa + sin2#sina(l — a)/2] 

— = do t (10) 

r [sin 29 cos a(l — a)/2 — sin a(sin 9 + a cos 2 9)] 

by recalling the relationship between Cartesian and circular coordinates. We 
note that when a = 1 (i.e in the circular-symmetric case) this equation reduces 
to 

* = d0^2 (11) 
r sin a 

and the solution is a line in the shape of a logarithmic (circular) spiral [3j . Solv- 
ing equation JHJ is a little more difficult because the presence of the coefficient 
a breaks the polar symmetry. We can start by studying two particular cases, 
i.e, when a = and a — tt/2. In the first case (a = 0) equation becomes 

dr cos 2 9 + a sin 2 9 

V = (a-l)sm9cos9 { } 

This equation can be solved, by rewriting it in the following form 
dr cos 9 asin# 

V = (a -1) sin 9 + (a-l)cos0 ( ' 



3 



and it gives 



where R is an integration constant. Equation i|14|) describes open lines, which 
tend to lines through the origin when a — ► 1, and that represent the analogous of 
radial-circular polarization for the geometry in subject, namely a radial-elliptical 
polarization. In Cartesian coordinates it reduces to the following relation be- 
tween x and y coordinates 

y = Kx a (15) 
where K is a constant value. In the other case (a = tt/2) equation (|10|) becomes 

£ = d9(a-l) , T 6C0S ° a „ (16) 

r v ' (1 + (a - 1) cos 2 9 v ' 

and it gives 

^ = fl (l + (a-l)cos^/ 2 ^ 

This is the polar form of ellipses, with (on considering a > 1) the major axis par- 
allel to y-axis and it represents the azimuthal polarization for elliptical geometry 
namely an azimuthal-elliptical polarization. When a = 1 the ellipses become 
circles (described by r = R relation), i.e. we find again the circularly-symmetric 
case. For the intermediate case, that is for a =/= 0, tt/2, the integration of the 
respective differential equation is less easy. In Figure we report the solution 
after having choosen a = 0.8 and the quotient a = B/A = 10/3. The field 
line that arises is a logarithmical (but elliptical) spiral. Spirally beams were 
introduced by Gori 0] for circular-symmetric case and here we found a natural 
extension to elliptically-symmetric case. In the following we analyse some prop- 
erties of this kind of beam. We have already seen which is the role of a angle: 
when it tends to zero the polarization tends to become radial-elliptical, the field 
being tangent to curves like in equation (|15J) whereas when it tends to oo the 
polarization tends to be azimuthal-elliptical, the field being tangent to ellipses 
described by equation Q17[l. Also we note that the beam described by equation 
l|5jl is anisotrope, i.e its tranverse intensity pattern depends by 9 angle. Indeed, 
for the intensity of transverse component holds 

I T (r,9,z) oc |E T | 2 =C 2 |/(r,z)| 2 [cosh 2 £-cos 2 (6> + a)] (18) 

in which the angular dependence is evident. It is simple to see that the maximum 
intensity directions are 

01 = 7r/2 - a ± rmr (19) 
(to = 0, 1, 2..) in which equation (|18fl becomes 

I T (r,6,z) oc |E T | 2 = C 2 |/(r,z)| 2 cosh 2 C (20) 

so we can find the value of angle a and, as a consequence, also the polariza- 
tion kind, by intensity measurement, contrary to the circular-symmetric case 
in which the intensity is dependent only from the radial coordinate r and does 
not give any information about polarization. In other words, in the circular- 
symmetric case all transverse intensity profiles associated to a particular polar- 
ization state (radial, azimuthal, spiral ones, etc.) are undistinguishable one to 
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another. In the elliptical-symmetric case, one can obtain all the parameters, i.e. 
the ratio B/A and the angle a, simply by the intensity transverse map measured 
across a generical tranverse plane. Indeed, with reference to equation l|18|) . by 
measuring the maximum intensity direction one evalues the angle a (equation 
(|19J) . whereas by measuring the angular- visibility ratio 

iMax ~ iMin 1 /ri -,N 

Vr = T = TJ7 (21) 

i-Max COSh £ 

one evalues the values of A and B (the indice r indicates that the radial coor- 
dinate is fixed). Also we note that the equation l|18|l never goes to zero, this 
because the term cosh 2 £ — cos 2 (9 + a) is different to zero V£ ^ (the case £ = 
represents linear polarization and here is not considered). It follows that it ever 
holds v < 1. It remains to see how one can obtain such a field. This is the 
purpose of the next section. 



3 Elliptically Symmetric Polarized Bessel-Gauss 
Beams 

It is well known that Bessel-Gauss beams (BG beams for short) are solution of 
scalar paraxial wave equation and were introduced for the first time from Gori 
et al. 0. In particular Gori introduced the linearly polarized BG beams of 
zero-order, here denoted for convenience as BGq. Subsequently many authors 
have been interested both to the generalization to vectorial 9 and higher-order 
BG beams, both to beam with unusual polarization features, as azimuthal and 
radial ones 6 . Here we recall, for reader convenience, the form of a BGq 

BG ° = (l+i/L) M (l + L/L) ] eM w^ + 1Z /L) ] eXp[ 2fc(^z/L) ] (22) 

where Jo is the zero-order Bessel function of the first kind, L is the Rayleigh 
length, w is the minimum spot size of the Gaussian beam, and (3 is the intensity 
of the tranverse wave vector; to extend the vectorial treatise to encompasse also a 
non-uniform elliptical polarization we start from a vectorial paraxial wave equa- 
tion. On starting from Maxwell's equations in vacuum, in absence of sources, 
we arrive to the following equation for the complete vector field E 

V x V x E = fc 2 E 

V • E = (23) 

On introducing a "slowly varying part" of E, say F, with 

E = Fexp(z£;z) (24) 

and by making use of the paraxial approximation, we obtain that F is governed 
from the following equations |1 0| 

V^Ft + 2ikd z F T = (25) 

V T ■ F T + ikF z = (26) 
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in which we have decomposed the field F in longitudinal and transverse com- 
ponent, i.e. F = Ft + F z z. Now, if we decompose also the field Fx into 
rectangular component, i.e. we let Ft = F x x + F y y we obtain two scalar parax- 
ial wave equations, each for single component, 

V 2 T F X + 2ikd z F x = (27) 

V%F y + 2ikd z F y = (28) 

equation (|27fl - l|28[l can be rewritten as LF X = and LF y = 0, after introduced 
a linear differential operator L = d x + d 2 + 2ikd z . Now, as the L operator 
commutes with the infinitesimal displacement operators [H] d x , d y and d z , it 
is possible to obtain new solutions from known solutions by mean of arbitrary 
differentiation. This observation will be crucial to obtain elliptically symmetric 
polarized BG beams. Indeed, the first thing to do is differentiating, respect to x 
and y variables, a BGq beam; in this way we obtain two new possible solutions 
for a scalar paraxial wave equation: d x BGo and d y BGo. If we now linearly 
combine them properly we find the ESPB. To this purpose it is sufficient to let 

F x = A[d x {BG ) cos a - d y {BG ) sin a] (29) 
Fj = B[d y (BG ) cos a + d x {BG ) sin a] (30) 

in which A and B are two arbitrary constant, with A ^ B, and a, as usual, is a 
constant angle. On utilizing the position in equations ©-0 wc finally obtain 

F T (r,z) = 

x 
+ 

As the paraxial propagation has to be satisfied, from equation (|31|1 follows that 
the function /(r, z) on the source plane [z — 0) can't be choosen in a arbitrary 
way. In particular, by comparing equation @ and equation <|31|) . we have 

f(r) = /(r,*)U=o = exp[-4][M(/?r) + *J o {0r)] (32) 

In Figures (0 we report the intensity pattern of a ESP Bessel-Gauss beam, 
evalued at different transverse planes. We see that the beam mantains, during 
its paraxial propagation, a zero in correspondence of z axis (a properties that is 
shared with the classical donut beams). In Figure we report the polarization 
pattern, in correspondence to the z = plane, in which is shown the behaviour 
of tranverse field for a — 0.8, and a ratio B/A = 10/3. On changing the sign 
of the linear combination in equation (|29|l as well as the ratio B/A one can 
obtain other interesting polarization pattern. In Figure (0J we plot intensity 
across the tranverse plane z = 0, for different values of angle a, to show the 
aforesaid angular dependence. Before closing the present section, we wish to say 
something about the longitudinal component of the electric field. To do that 
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■exp[ 



{1 + iz/L) 2 w 2 (l + iz/L)' 

, if3 2 z „, T , (3r 

cxph 2fc(i + *W [/3Jl( IT^7Z ) + 

2r fir 

~^^o(z n -)][sinh£cos(# + a)x + cosh£ sin(6' + a)y](31) 

w z 1 + iz/L 



we will refer to the generical form of EPSB as that in equation JSJ and we shall 
use equation (|26() to obtain the longitudinal component in paraxial regime: 

F z = ^Vr • F T (33) 
k 

After straightforward calculations we find 

F z = (d r f(r,z) + ^^)smh^cosa + exp(-md r f(r,z)sm 2 e+ (34) 
r 

+ ^^-cos 2 e)cosa+-(d r f(r,z) - Zillil) s ina sin 201 
r 2 r 

which represents the full form of the longitutinal component in paraxial approx- 
imation. Here we see that, contrary to the circular-symmetric case for which an 
azimuthally polarized field does not give rise to a component along the z-axis, 
when we have an azimuthal-elliptical polarization (for which a = 7r/2), F z is 
not null and holds 

F z = exp(-a^(d,/(r,z) - ^M)sin20] (35) 

and it goes to zero as exp(— £). This is obvious if we recall that when £ — > 00 
the geometry polarization tends to become circular because for high value of £, 
sinh£ « cosh£ « exp(£) and the ellipses in equation (II Tf) tend to become cir- 
cles. The fact that the longitudinal component changes by changing coordinate 
system is well-known and it can explained by mean a geometrical point of view. 
Indeed, equation (|33|l states that F z is a divergence of a vectorial field. The 
divergence Vt ■ Ft counts the algebrical sum of sinks (i.e. where a field line 
of Ft terminates) and sources (i.e. where a field line of Ft originates) into a 
surface element, across a (x, y) plane. Since the divergence is a pseudoscalar, 
namely we are dealing with a density, it is not invariant under a generical space 
deformation. If, for example, the space is squeezed in a way in which the areas 
are reduced, the density (i.e. F z ) increases. In practice, a distortion of space can 
be thought as a transformation of coordinates and the fact that the divergence 
operation is dependent by the coordinate system in which we work means that 
we are dealing with metrical, and not topological, features of the field [TT] . 



4 Conclusion 

We have introduced unconventionally polarized optical beams endowed with 
elliptically-symmetric features in their tranverse polarization pattern. A beam 
of this kind mantains its polarization state under free-space paraxial propaga- 
tion. In particular we have derived a natural extension to elliptically-symmetric 
case of spirally-polarized beams introduced by Gori 0j for circular symmetric 
case. Some properties of this kind of field have been analysed. In particular, we 
have pointed out that, by using an elliptically-symmetric polarization, the tran- 
verse intensity profile gives directly information about the polarization state, 
the angular behaviour of the intensity depending from it. This implies that it 
is possible to measure the field polarization state only by an acquisition of the 
intensity pattern across a generical tranverse plane. This is in contrast with the 
analogous situation in circular-symmetric case, in which the intensity pattern is 
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indipendent from the particular polarization state. Indeed, a field polarized in 
a radially way produces a tranverse intensity profile identical to that produced 
by the azimuthally or spirally polarized ones, or any other combination of them 
which mantains the circular symmetry. Due to its important role in optics mi- 
croscopy, we have reported also the form of the longitudinal component for the 
electric field, in paraxial approximation. 
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Figure 1: elliptical logarithmic spiral. The parameters are choosen as follows: 
B/A = 10/3, a = 0.8. 
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Figure 2: Tranverse intensity of a ESP Bessel-Gauss beam with non- uniform 
elliptical polarization at different (x,y) plane, a) z — 0; b)z — 0.1L;3)z = 
0.2L;4)z = 0.3L, where L is Rayleigh length. The intensity is normalized at 
one to each plane. With reference to equation we have f3 = 8/w, a = 0.8, 
B/A = 2. The axises are normalized to w value. 
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Figure 3: Polarization pattern of a ES?lBcssel-Gauss beam across source plane 
(z = 0). The parameters are choosen as follows: B/A = 10/3, (3 = 8/w and 
a)a = 0, b)a = ir/2, c)a = 0.8. The axises are normalized to w value. 




Figure 4: Intensity map of a ESP BesseJ-*pauss beam across source plane (z = 0) 
for different values of polarization angle a: a) a = 0, b) a = ir/6 c) a = ir/A 
d)a = it/2. The parameters are choosen as follows: B/A — 10/3,/? = 8/w. The 
axises are normalized to w value. 



